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Abstract 

Endometriosis is one of the most frequent benign gynecological disorders. Numerous studies have shown an association 
between GSTM! and/or GSTT1 polymorphisms and endometriosis susceptibility. However, these associations remain 
inconclusive. To derive a more precise estimation, we conducted a comprehensive search to identify all existing studies and 
then performed a meta-analysis. Electronic literature searches of the PubMed, Chinese Biomedical, and China National 
Knowledge Infrastructure databases were performed up to December 2013. GSTM1-, GSTT1-, and dual-null genotypes were 
analyzed independently, and pooled odds ratios (ORs) with 95% confidence intervals (95% CIs) were calculated by 
comparing the null genotype with other genotypes using the random-effects or fixed-effects model. Twenty-five and 16 
independent studies on GSTM1 and GSTT1 polymorphisms, respectively, and five GSTM1-GSTT1 interaction analyses were 
identified and included in this meta-analysis. Both GSTM1- and GS7T7-null genotypes increased risk of endometriosis 
(OR=1.54, 95% CI: 1.30-1.83, P<0.001; OR=1.41, 95% CI: 1.10-1.82, P=0.007; respectively). Moreover, we found a 
significant positive association between the dual null genotype GSTM1-GSTT1 and endometriosis susceptibility (OR=1.33, 
95% CI: 1.03-1.72, P = 0.027). This meta-analysis provides evidence that null genotypes of GSTM1 and/or GSTT1 contribute to 
risk of endometriosis. Further investigations are required to confirm these findings. 
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Introduction 

Endometriosis, a benign gynecological disease, is characterized 
by the presence of endometrial glands and stroma at extrauterine 
sites. Approximately 6-10% of women of reproductive age suffer 
from this condition [1], [2], which negatively impacts quality of life 
by causing pelvic pain, heavy menstrual flow, dysmenorrhea, and 
infertility [3], [4], making endometriosis a major public health 
threat. 

The etiology and pathogenesis of endometriosis remain unclear. 
Endometriosis is commonly considered as a complex trait caused 
by the interaction between genetic and environmental factors [5] . 
Both genetic polymorphisms and environmental factors are 
considered risk factors for endometriosis [6], [7]. Even though 
no simple Mendelian inheritance was confirmed, geneticfactors 
increased the risk of endometriosis by 6% for near relatives [8]. 

Polymorphisms in the glutathione S-transferase (GST) system 
have long been recognized as a risk factor for endometriosis and 
have been extensively explored. Human GSTs are a multigene 
family of phase II metabolizing enzymes that are crucial in 
detoxification of xenobiotics such as carcinogens, environmental 
toxins and drugs [9]. Human cytosolic GSTs have been 
subdivided into eight distinct classes: alpha (GSTA), mu (GSTM), 
pi (GSTP), theta (GSTT), kappa (GSTK), zeta (GSTZ), omega 
(GSTO), and sigma (GSTS) [10]. Among the genes in the GST 



superfamily, GSTM1 (lpl3.3, MIM 138350) [11] and GSTT1 
(22qll.2, MIM 134660) [12] are the most extensively studied 
owing to their critical role in detoxification and the high-frequency 
of allelic variants. GSTM1 and GSTT1 polymorphisms have been 
suggested to be associated with endometriosis by many epidemi- 
ological studies [13]-[24]. However, findings on the direction of 
the association remain equivocal. 

A previous meta-analysis conducted by Sun-Wei Guo [25] in 
2005 indicated a significant association between the GSTT1 -null 
genotype and endometriosis, however no such association was 
found between the GSTM1- or GSTM 1 -GSTT 1 -null genotypes 
and endometriosis risk. Since then, 12 relevant studies [5], [8], 
[19]-[24], [26]-[29] have further examined the associations 
between the two polymorphisms and endometriosis risk. We 
aimed to confirm the potential associations by conducting an 
updated meta-analysis, to provide insight into the pathophysiology 
of endometriosis. 

Methods 

Identification and eligibility of studies 

Studies were identified by searching the PubMed, CBM 
(Chinese Biomedical), and CNKI (China National Knowledge 
Infrastructure) databases for relevant reports published prior to 
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December 2013, using the key words "GSTM1", "GSTT1", 
"polymorphisms", "endometriosis", and combined phrases. Ad- 
ditional literature was collected from cross-references within both 
original and review articles. No restrictions on language, 
population, or sample size were set in this meta-analysis. Studies 
were required to comply with the following inclusion criteria: (1) 
original case-control or cohort studies; (2) studies investigated the 
association of GSTM1 or GSTT1 polymorphism with risk of 
endometriosis; (3) sufficient information to calculate odds ratios 
(ORs) with 95% confidence intervals (CIs); and (4) Chinese articles 
were published in Chinese core periodicals. Exclusion criteria 
were: (1) not case-control or cohort studies evaluating the 
association of the GSTM1 or GSTT1 polymorphism with 
endometriosis; (2) case reports, letters, reviews, editorials, or 
correspondence articles; (3) studies based on incomplete raw data; 
and (4) studies that contained overlapping data. 

Data extraction 

The data from eligible studies was independently checked and 
extracted according to the pre-specified selection criteria, and the 
discrepancies were resolved by discussion and agreement between 
all investigators. The following information was collected from 
each included study: name of the first author, publication year, 
study location, ethnicity, source of controls, sample-size of cases 
and controls, and genotype frequency in cases and controls. 
Different ethnicity descents were categorized as European and 
Asian. According to source of controls, all included studies were 
defined as population-based (PB) and hospital-based (HB). 

Quality score assessment 

Two investigators independendy assessed the quality of included 
studies using the Newcasde-Ottawa Scale (NOS) [30] . The NOS 
ranges from zero to nine stars. Studies with a score of seven stars 
or greater were considered to be of high quality. Discrepancies 
were resolved as described above. 

Statistical analysis 

Meta-analyses were performed for polymorphisms investigated 
in at least three studies [31], [32]. The strength of the associations 
between GSTM1 and GSTT1 polymorphisms and endometriosis 
risk were measured by ORs and respective 95% CIs, which were 
calculated by comparing the null genotype with other genotypes. 
The significance of the pooled ORs was determined by the Z-test 
(P<0.05 was considered significant). Subgroup analyses were 
performed by ethnicity and source of controls. 

Heterogeneity among studies in terms of degree of association 
was evaluated using % 2 tests. The I 2 statistic was used to evaluate 
variations due to heterogeneity rather than chance. P<0.10 or 
I 2 >50% indicated the presence of between-study heterogeneity, 
and the random-effects model (DerSimonian-Laird method) [33] 
was used to calculate the pooled ORs; otherwise, the fixed- effects 
model (Mantel-Haenszel method) [34] was selected. 

Sensitivity analyses were performed to evaluate the stability of 
the results of the meta-analysis. The influence of individual studies 
was evaluated by estimating the pooled ORs in the absence of 
each study [35]. Potential publication bias was investigated by 
visual inspection of Begg's funnel plots. We also used the Begg's 
[36] and Egger's tests [37] to evaluate any possible publication bias 
(P<0.05 was treated as significant publication bias). All statistical 
analyses were performed using Stata statistical software, version 
12.0 (Stata Corp., College Station, TX, USA). 



Results 

Study selection and characteristics 

This meta-analysis is guided by the PRISMA statement 
(Checklist SI). We initially identified 85 results relevant to the 
search terms in the selected databases. After reading the titles and 
abstracts, only 28 potentially eligible articles were identified for 
further detailed evaluation. Of these, two articles [38], [25] were 
excluded because one article was a review and the other was a 
meta-analysis. After further screening, two further articles [39], 
[40] were excluded for overlapping populations. Two ethnic 
populations (Hans and Uyghurs) were studied in the article by 
Ding et al [18], so we considered it as two studies in this analysis. 
In total, we included 24 articles [5], [8], [13]-[24], [26]-[29], 
[41]— [46], including 25 independent case-control studies (Fig- 
ure 1). Eighteen, four, and two of the articles were written in 
English, Chinese, and Russian, respectively. Twenty-five and 16 
studies focused on GSTM1 and GSTT1 polymorphisms, respec- 
tively, and five were on GSTM1-GSTT1 interaction analysis. The 
characteristics of the included studies are described in Table 1. 

Quality assessment results 

The scores of included studies were 5 to 9 (Table 1). 

Meta-analysis results 

The association between GSTM1 polymorphism and endome- 
triosis was investigated in 25 studies, which included a total of 
3330 cases and 3959 controls. The heterogeneity was significant, 
so the random-effects model was selected. The result showed that 
the null genotype of GSTM1 was associated with an increased 
endometriosis risk (OR= 1.54, 95% CI: 1.30-1.83, P<0.001). A 
forest plot is shown in Figure 2. Furthermore, we included 22 
studies of high quality to validate the association, and the result 
again showed a strong association between this polymorphism and 
endometriosis risk (OR= 1.45, 95% CI: 1.23-1.71, P<0.001). 
Subgroup analysis by ethnicity was performed, and an increased 
risk of endometriosis was observed in Europeans and Asians 
(OR= 1.58, 95% CI: 1.19-2.09, P = 0.002; OR= 1.52, 95% CI: 
1.21—1.91, P<0.001; respectively). In stratified analysis by source 
of controls, the results showed that source of controls did not affect 
the pooled results and a significantly increased risk of endometri- 
osis was detected both in PB and HB studies (OR= 1.52, 95% CI: 
1.08-2.16, P = 0.005; OR= 1.55, 95% CI: 1.26-1.91, P<0.001; 
respectively). 

Sixteen independent studies, with a total of 2263 cases and 2380 
controls, were included in the meta-analysis of GSTT1 polymor- 
phism. We found significant heterogeneity between studies, so the 
random-effects model was used to pool the results. The results 
indicated a positive association between the null genotype of 
GSTT1 with endometriosis risk (OR= 1.41, 95% CI: 1.10-1.82, 
P = 0.007). A forest plot is shown in Figure 3. Stratified analysis 
by ethnicity showed a significant association in Asians (OR= 1.53, 
95% CI: 1.14-2.06, P = 0.005), but not in Europeans (OR= 1.21, 
95% CI: 0.74-1.99, P = 0.45). 

For GSTM1-GSTT1 interaction analysis, five independent 
studies, including 1008 cases and 1 164 controls, were subjected to 
meta-analysis. The heterogeneity was not significant, so we 
selected the fixed-effects model. The results showed that the dual 
null genotype of GSTM1-GSTT1 was associated with an 
increased endometriosis risk (OR= 1.33, 95% CI: 1.03-1.72, 
P = 0.027). A forest plot is shown in Figure 4. The results of this 
meta-analysis are summarized in Table 2. 
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85 potentially relevant studies identified from PubMed, 
CBM and CNKI databases after initial search 



30 articles were excluded for they were identical articles. 
27 articles were excluded for not case-control studies 



28 articles concerning the current topic for full-text view 



1 article was excluded for it was review; 

1 article was excluded for it was a meta-analysis 



26 articles relevant to GSTMlovGSTTl genotypes and 
endometriosis for detailed evaluation 



2 articles were excluded for they were reduplicative studies 



Finally, 24 articles included in our meta-analysis 



Figure 1. Flow diagram of the study selection process. CBM: Chinese Biomedical; CNKI: China National Knowledge Infrastructure. 
doi:1 0.1 371 /journal.pone.01 06761 .g001 



Sensitivity analysis 

Sensitivity analyses were performed after the sequential removal 
of each eligible study to assess the influence of each individual 
study on the pooled ORs. In the analysis of the GSTM1 
polymorphism, the pooled ORs were not qualitatively changed 
when any single study was omitted, suggesting that no single study 
exhibited excessive influence, and that the results are reliable 
(Figure 5). Other results were also relatively stable. 

Publication Bias 

We conducted Begg's and Egger's tests to evaluate potential 
publication bias. There was no statistical evidence of publication 
bias regarding analysis of the GSTM1 polymorphism (GSTM1 (— ) 
vs. GSTM1 (+): Begg, P = 0. 129 and Egger, P = 0.079), and Begg's 
funnel plots suggested no substantial asymmetry (Figure 6). 
There was no publication bias for other results (GSTT1 (— ) vs. 
GSTT1 (+): Begg, P = 0.344 and Egger, P = 0.207, and the shape 
of the funnel plot also did not reveal any evidence of obvious 
asymmetry (Figure 7); GSTM1-GSTT1 (— ) vs. other genotypes: 
Begg, P = 0.462 and Egger, P = 0.613). 

Discussion 

As many GST genes are polymorphic, whether particular allelic 
variants in GST genes are correlated with altered risk of some 
kinds of diseases has provoked great interest [47]. The null 
genotypes of GSTM1 and GSTT1, two of the most widely-studied 
polymorphisms, are characterized by homozygous deletions of the 
respective genes [39]. The first study considering null genotype of 
GSTM1 as risk factors was conducted by Baranov et al. [13] in 



1996. A preponderance of GS7T7-null subjects among endome- 
triosis patients was detected in 1999, although it was not 
statistically significant [14]. Since then, many studies have 
investigated the associations between null genotypes of GSTM1 
and GSTT1 and endometriosis susceptibility. However, the results 
are inconsistent and conflict, which compelled us to pay attention 
to the two polymorphisms at a meta-analytical level. 

The current study, including 25 case-control studies with 3330 
cases and 3959 controls, is the most comprehensive meta-analysis 
of the association of GSTM1- and GSTT1 -null genotypes with 
endometriosis risk, which allows us to expand the discussion of 
possible implications and interpretations of the findings. Our 
meta-analysis suggested that there were significant associations 
between the null genotypes of GSTM1 and GSTT1 and 
endometriosis. Moreover, the combined GSTMl-GSTTl-ixaR 
genotype also showed a positive association with endometriosis 
susceptibility. 

A novel finding of the present study was the significant positive 
association between the null genotype of GSTM1 and endome- 
triosis risk. The large number of studies and subjects included in 
this meta-analysis were substantial enough to resolve the issue of 
conflicting results obtained in individual studies, which was 
primarily caused by their small sample size. Furthermore, the 
results in subgroup analysis by ethnicity also indicated an 
association in both Europeans and Asians suggesting this 
association is reliable. 

This meta-analysis showed a moderate positive association 
between the null genotype of GSTT1 and endometriosis risk, 
which is consistent with the result of a previous meta-analysis [25]. 
In the analysis stratified by ethnicity, a significant association was 
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ID 
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Figure 2. Forest plot of pooled OR with 95% CI for association between the null genotype of GSTM1 and endometriosis risk. 

doi:1 0.1 371 /journal.pone.01 06761 .g002 



found in Asians, but no such association was detected among 
Europeans. There are several possible reasons for such a 
difference. First, the frequencies of the risk-associated homozygous 
null genotype vary between different races. The frequency of the 
GST77-null genotype is nearly 50% in the Chinese population 
[48], [49] 14.5-20.1% in Indians [50]-[52], and 11.0-37.9% in 
Europeans [53], [54]. Thus, the GSTT1 polymorphism may exert 
varying effects in different populations. Second, the different 
results could also be explained by study design or sample size. 
Other confounding factors, such as age and lifestyle may also be 
considered. 

The GSTM1-GSTT1 interaction analysis indicated that women 
with double-null genotype had significantly increased endometri- 
osis risk compared with those with other genotypes. If genetic 
susceptibility to endometriosis is, at least in part, mediated through 
polymorphisms of genes that encode enzymes responsible for 
detoxification, it is possible that the combination of GSTM1- and 
GST77-null genotypes may be more discriminating as a risk factor 
for endometriosis than a single null genotype. 

The observation that the null genotypes oiGSTMl, GSTT1, 
and GSTM1-GSTT1 increased the risk of developing endometri- 
osis is biologically plausible. Environmental contaminants, such as 
polychlorinated dibenzo-p-dioxins, polychlorinated biphenyls, and 



polycyclic aromatic hydrocarbon, have been suggested to promote 
the occurrence and development of endometriosis by interfering 
with the estrogen signaling pathway and their immunosuppressive 
effects [55] . The GSTs plays a critical role in the detoxification of 
a broad range of environmental contaminants. As critical phase II 
metabolic enzymes, GSTs catalyze reactions between glutathione 
and all kinds of potentially lipophilic compounds, causing 
neutralization of the carcinogens, products of oxidative stress 
and toxic compounds [56]. Previous studies suggested that 
homozygous null deletions in GSTM1 and GSTT1 cause a 
complete loss of the activity of their encoded enzymes [57], [58]. 
The GSTM1 and GSTT1 deletions are detected in 42-60% and 
13-26% of Caucasians, respectively [59]. Lack of GST enzyme 
activity resulting from the null genotypes may affect the 
detoxification of environmental toxins, and thus contribute to 
the pathogenesis of endometriosis. 

Heterogeneity is an important issue in meta-analysis. Although 
we minimized the likelihood by performing a careful search for 
published studies, using the explicit criteria for study inclusion, 
statistically significant heterogeneity still existed in most compar- 
isons. There are several explanations for the significant between- 
study heterogeneity, such as different study populations, genetic 
factors, and environmental factors. In particular, environmental 
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Figure 3. Forest plot of pooled OR with 95% CI for association between the null genotype of GSTT1 and endometriosis risk. 
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Figure 4. Forest plot of pooled OR with 95% CI for association between the dual null genotype of GSTM1-GSTT1 and endometriosis 
risk. 
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Figure 5. Sensitivity analysis for the meta-analysis regarding the association between GSTM1 polymorphism and endometriosis 
risk. 
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contaminant exposure, a risk factor for endometriosis, is an 
important factor contributing to heterogeneity. The status of 
environmental contamination varies between countries, so the 
endometriosis incidence varies between populations. Publication 
bias is another important issue which should also be accounted for 
in meta-analyses. After evaluating the publication bias using 
Begg's funnel plots we did not detect a publication bias, indicating 
the strength of the results. 

Sun-Wei Guo [25] also evaluated the association between 
GSTM1 and GSTT1 polymorphisms and endometriosis risk by 
performing a meta-analysis including 14 studies with 1539 cases 



Begg's funnel plot with pseudo 95% confidence limits 
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Figure 6. Begg's funnel plot for publication bias in selection of 
studies regarding the GSTM1 polymorphism. 
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and 1805 controls. That study suggested that the endometriosis 
risk associated with the null genotype of GSTT1 was 29% higher 
than other genotypes, but it failed to find positive associations 
between the null genotype oiGSTMl or GSTM1-GSTT1 and the 
risk of endometriosis. There were some differences between that 
study and ours. First, our meta-analysis provided more compre- 
hensive information on the relationships between the two 
polymorphisms and endometriosis by extracting data from more 
studies with more total cases and controls. Second, some issues 

Begg's funnel plot with pseudo 95% confidence limits 
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Figure 7. Begg's funnel plot for publication bias in selection of 
studies regarding the GSTT1 polymorphism. 
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that may affect the results of meta-analysis, such as publication 
bias, sensitivity analysis, and quality assessment of the included 
studies, were addressed in our study. Third, the current study also 
showed distinct findings, with the added statistical power. 

This meta-analysis had several limitations that should be taken 
into account when considering its contributions. First, heteroge- 
neity among studies existed in some comparisons of polymor- 
phisms. Second, the results of our meta-analysis were applicable to 
only two ethnic groups, Europeans and Asians, as there were no 
relevant studies including data from African ethnic groups. Hence, 
to conduct a more precise analysis of the association between 
GSTM1- and GSTTl-null genotypes and endometriosis risk, 
additional studies with larger sample sizes and involving different 
ethnicities (especially African) are needed. Third, gene-environ- 
ment interactions were not evaluated in this meta-analysis. 
Subgroup analyses based on environmental exposures were not 
performed because of insufficient data on such associations in all 
included studies. 
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